The influences of long range Coulomb interaction on the electronic Mach-Zehnder 

interferometer of quantum Hall edge states 



Hyun C. Lee 

Department of Physics and Basic Science Research Institute, Sogang University, Seoul, Korea 

(Dated: April 23, 2009) 

The influences of long range Coulomb interaction (LRCI) on Mach-Zehnder interferometer (MZI) 
constructed on quantum Hall edge states is studied employing bosonization method. The interaction 
of interchannel zero-modes is shown to give rise to a characteristic energy scale which is of the order 
of the period of experimentally observed lobe pattern of visibility. The nonmonotonic behavior of 
visibility as found Chalker et al. is understood analytically using asymptotic analysis. 
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Introduction- An electronic MZI is a type of matter 
wave interferometer that has been reahzed with (integer) 
QH edge systems jii^ii^ Among many interesting phenom- 
ena observed in electronic MZI we focus on the puzzling 
lobe pattern of the interference visibility of differential 
conductance a{V) {V is a bias voltage) 

This lobe pattern is hard to understand in the frame- 
work of non-interacting electrons so it is generally be- 
lieved to be due to many-body interaction. A few the- 
oretical proposals have been made to explain this lobe 
pattern: the introduction of additional edge modea^i^, 
the decoherence and finite temperature effect^, the shot- 
noise efFects&. The visibility of MZI of fractional QH 
edges which does not include LRCI has been also studied 
recently.^ 

In this Brief Report we generalize the approach of 
Ref. [3| in two ways: (1) LRCI of the zero-modes (see 
below) between two channels which comprise MZI is in- 
cluded (2) both the integer QH states and the Laughlin 
fractional QH states at filling fraction v = l/(2n-|-l), n = 
1,2,... are considered. In the fractional case we have to 
take both the fractional quasiparticle and the electron 
tunneling into account at point contacts. 

A possible relevance of interchannel LRCI may be ar- 
gued as follows. The length scale of MZI in Ref.01 is 
about R 5/im. Taking the typical dielectric constant 
of QH devices to be e ~ 10, the associated interchan- 
nel Coulomb energy scale Ec is / Re ~ 5 x 10~^/zeV, 
which is of the same order of magnitude as the observed 
period of lobe pattern of visibility (see Fig. 3 and Fig. 4 
in Ref.jl]). This energy scale is well expected to depend 
on the geometry of MZI as well as the applied magnetic 
field. The investigation of this interchannel interaction of 
zero-modes (defined below) requires a very careful treate- 
ment of the so-called Klein factors of bosonization for- 
mula, which constitutes the most significant part of this 
report. 

The results of Ref. [7|, in particular, the nonmonotonic 
behavior of visibility as a function of bias voltage (Fig. 6 
of Ref. ^ ) , depend crucially on the asymptotic behavior 
of one-electron Green's function in the presence of many- 
body interactions. We analyze the asymptotic behavior 
of Green's function employing the method of asymptotic 



analysis^*' , thus providing more analytic understanding 
of the results. It turns out that the momentum depen- 
dence of interaction matrix element plays an important 
role as we will discuss. The main results of this report are 
(1) the exact time-dependent Klein factor which give rise 
to interchannel Coulomb energy scale Eq.([T6]);(2) the 
tunneling current which incorporates the energy scale E^^ 
Eqs. (l23l24p :f3) the analytic form of the visibility Eg. ([55)1 . 

Model of MZI of QH edges - The Hamiltonian which 
applies to both the integer and fractional Laughlin QH 
edge state and acts within each channel isiiiiS 

Hintra = ^ — dx [pi{x)]^ 
i=l,2 

+ i / dxdyV{x-y) ^ p^{x)p^{y), (1) 

i=l,2 

where i = 1, 2 is the channel index of MZI (see Fig.l of 
Ref. 7] for a schematic view of MZI) . pi{x) — + ^ dx4>i 

is the density operator of edge i. Ni is the number op- 
erator of edge i whose momentum is zero, hence it is of- 
ten referred to as a zero-mode operator in the context of 
bosonization>i2, i/ is the filling fraction of QH system and 
Lg is the system size. 4>i is the boson operator which de- 
scribes the collective harmonic modes with nonzero mo- 
mentum. Vi is the velocity of collective modes of chan- 
nel i. Two channels will be assumed to be identical, so 
that vi = V2 = vq. V{x) = ^^^'i^^i is the LRCI acting 
within each channel, a is a short distance cutoff. The 
fundamental commutation relation of density operators 
is" [pi{x),pj{y)] ^ihSij^dxS{x-y). 

There are two interactions which couple two channels 
i ~ 1,2. One is the tunneling interaction at two quantum 
point contacts (a and b) which play the role of beam 
splitters of an optical interferometer. 

Ht = ia*l(0)*2(0) + 4*l(x - /i)*2(x = h) + h.c, (2) 

where the operator 5'i(x) can be either electron or quasi- 
particle operator depending on the character of the point 
contacts. The bosonized expression of these operators are 
given in Eqs. (|7llip . The tunneling amplitudes depend on 
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the enclosed flux $ via tJl = l^atfe |e'*/*° , where $o is 
the flux quantum J 

The other interaction which can couple two channels is 
LRCI between two channels. For simplicity we consider 
the interaction between zero-modes only. This is because 
the intcrchannel LRCI becomes singular logarithmically 
in the zero momentum limit (the singularity is cut by the 
finite system size) . The Hamiltonian for the interchannel 
LRCI is taken to be 



Hc = EcNiN2. 



(3) 



Extracting the zero-mode parts from Eq.([T]) and Eq.Q 
we can define a Hamiltonian for zero-modes only. 

Hzcro = ^ ^ + Eintr. + ^-^1^2 • (4) 



-Eintra IS an energy scale from the intra channel Coulomb 
intraction (the second term of Eq.((T|)). The chemical 
potential of each channel is determined by the average 
number of electrons in each channel {Ni)^ Let us de- 
fine an operators 5Ni which describes the fluctuation of 
electron number around the average value. 



8N, = N,-{N,), {6N,) =0. 



(5) 



The Hamiltonian for the fluctuations 5Ni can be obtained 
by inserting Eq.® into Eq.ljl]) and discarding constant 
terms. 



Ho = ^ ^Ji^5N,+Ec5mN2, 



(6) 



4=1.2 



where a contribution from i?intia Yli is neglected since 
the relative fluctuations between two channels will play 
a more important role in interference. As will be shown 
below, with the form of Eq.®, the time evolution of 
Klein factors can be determined exactly. 

Bosonization of electron and quasiparticle operators- 
The bosonized expression for the electron operator in the 
edge state of Laughlin quantum Hall liquid at filling frac- 
tion ly is given byiiii^ii^ {i — 1,2) 



*o,*(a;) 



1 



where i^o,i is the Klein factor which implements the fermi 
statistics of electron operators of different species. It 
satisfies the following relations. 



-P -A- ■ F F''' — 1 



{Fe,„ Fe,,} = 0, {Fe,„ F^} = 0, I ^ J. (8) 

Using the bosonization formula of Ref.jl^ (with appro- 
priate changes of notations) Eq. ([8]) can be explictly real- 
ized as follows: 



Fo 1 = e*''(^i+^2)/2i^ e^^^i/'', 



(9) 



where 0i 2 are the operators with zero momentum which 
are dual to Ni in the following sense. 

[e„N,] = +wS,,. (10) 

The bosonized expression of the quasiparticle operator 
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(2^)1/2 



(11) 



We could not find the explicit expression for the Klein 
factors of quasiparticle operator in literature. Based on 
the idea that the quasiparticle behaves like a fraction of 
an electron, we can take the (l/z^)-th root (recall 1/;^ is 
an integer) of Eq.®, thus leading to: 

Fq 1 = gi^(JVi+JV2)/2 g-ifli ^ 

p^^^^e-'-iNl+N2)/2^-^9,^ (^2) 

The validity of Eq. p2|) can be confirmed by the fact that 
Fq.i and Fq^2 satisfy the following commutation relations 
of fractional statistics. 



[7V„Fq,,] =-Z.Fqj%, Fq,lF, 



q.2 



e-""i^q,2Fq4- (13) 



Time evolution of Klein factors- The time evolution of 
Klein factors under the action of the Hamiltonian Eq. ^ 
can be determined exactly. We use the following opera- 
tor identityJ^ Let A,B,D be some operators satisfying 
[A, B] DB, [A, D] = [B, D] = 0. Then for arbitrary 
function f{A) of the operator A, we have 



f{A)B^Bf{A + D). (14) 
-Ail - E,5N2 (15) 



Idenftifying 

A^Ho, B^F.^u D 
we find (with f{A) = e*^*) 

Fe,i(i) = Fe,i(i = 0) e-Vlt-«Bc5Ar2t^ (^g) 

The dependence of Fe^i(t) on SN2 gives rise to additional 
time dependence for correlation functions, which is not 
present in Ref.F^. This time dependence will make the 
visibility exhibit features around the energy scale Ec [see 
Eq.dSSD]. Similarly, 



FeMt) = Pe,i{t = 0) e-*^^*-*-^=<='5^i*. (17) 
As for quasiparticle operators, we have 

Fq,l(i) = Fq,l(i = ^-i'^f^lt^i'^E^SN^t ^ 

Fq.2{t) = Fq,2(i = 0) e-^'''^^*-^''^-''^!*. (18) 

Tunneling Currents- The tunneling current operator is 
given by (electron tunneling is assumed) 



dN2 
dt 



total 



N2 



(ia*l(0)*2(0) + <b*l(/l)*2(/2) - h.c.)l9) 
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where i?totai is the sum of Eas. (|ll2l3p . A straightforward 
time-dependent pertm'bation in f, yields^ 



and the flux dependent part is 



/EE(/(i = 0)) 



- / dt(Y,Xk{t) +h.c. 



Xi{t) = |t,|2([vl/t(0,0)*2(0,0),vl.^(0,i)*i(0,i)]), 

X^(t) - tatl{[^\{0, 0)^2(0, 0), ^\{l2,t)-^l{llM, 

X4t) - t:u([^\{h,0)^2{l2,0), *k0, i)*i(0, t)]), 



(20) 



where the expectation value should be taken with respect 
to the eigenstates of the Hamiltonian Hintra + He [Eqs.dU 
[5]. For the quasiparticles, the current / must be multi- 
plied by v. The visibility of interference originates from 
the correlation functions X3{t) and X4{t). 

Each correlation function Xi{t) can be factorized into 
the product of the nonzero-mode (also called oscillators) 
contribution and the zero-mode contribution. Further- 
more the nonzero-mode contribution is a product of that 
of channel 1 and channel 2, since the Hamiltonian of the 
nonzero-modes of channel 1 and 2 commutes each other. 
The nonzero-mode contribution is identical with that of 
Ref.Q, and the details are omitted. The zero-mode con- 
tribution, however, cannot be expressed as a product of 
contribution from each channel, since the zero-modes of 
two channel are coupled by the Hamiltonian Eq.([3]). 

In spite of this noncommutativity, the zero-mode con- 
tribution can be computed exactly owing to the results 
Eqs. (ll6|18|) . For example, the zero- mode contribution of 
Xi{t) in the case of the quasiparticle tunneling is given 
by 



{Fl{0)F^2{0)FUt)F^iit)) = e 



(21) 



To obtain Eq.(ini) we insert the rcsuhs Eq. (|16ll8p for the 
time dependent Klein factors, then employ the commu- 
tation relation between Klein factor and the zero-mode 
operator 6Ni, and use {SNi) — 0. 

Using complex conjugation properties of Xi it can be 
shown that (/q and defined in an obvious way) 



dt 



Yoit) + 2ReY^{t) . (22) 



The visibility is basically determined by the flux depe- 
dendent conductance f7$ = dIq,/dV {eV = /X2 — /ii). For 
electron tunneling, the flux independent part is 



i\ta\ 



\tb 



2\ i(M2-Ml)t 



(27ra)2 



h.c. 



(23) 



b gi(/i2-A'i)t ^2'iTi({N2)l2-{Ni)li)/vLs 



" (27ra)2 



(24) 



ge{x,t) is the nonzero- mode contribution for the one- 
electron Green's function, which has been obtained in 
Ref.^. At T = 0, 

5e(a;, t) = exp - Ce(x, t) - ise{x, t) , 

Ce(x,t) = -/ —er'"^(l-cos(qx + uj„t)], 
Jo 1 ^ ' 



Se{x,t) = — [ —e '^'^ sm{qx + uj„t), 
V Jo q 

where the frequency of nonzero-modes is given by 

i^q^q[vo + Vcm — , Vc = 
qa 



27rfe' 



(25) 



(26) 



where ga < 1 is assumed and C ~ 1.13 is a numerical con- 
stant. The logarithmic factor of Eg. ipS)) comes from the 
intrachannel Coulomb interaction matrix element. The 
major difference from those of Ref. [?| is the presence of 
the factor e^"^^"*" which originates from the interchannel 
LRCI. For quasiparticle tunneling, the following modifi- 
cations are to be made 



Ec^iy^Ec, (Ai2 - Ml) ^^(^2 - 

Ceix, t) V^Ce{x, t), S^ix, t) -> v'^Se{x, t) . 



(27) 



We note that the energy scale enters in such a way 
it is not a mere additive renormalization of chemical po- 
tential. Also, in higher 2n-th order expansion in tunnel- 
ing amplitude iai^fc: a factor e^"'-^"'^ will emerge since 
this factor is generated by commutation of Klein factors. 
Recall that the Klein factors accompany the tunneling 
amplitude in Eq.([21). Thus, one can expect some features 
of visibility which is due to geix, t) will appear being cen- 
tered at energies nEc with decreasing magnitude. In the 
second order expansion of ta,b, the visibility (or equiv- 
alently flux dependent conductance (T$) is found to ex- 
hibit nonmonotonic behaviori^ Let us try to understand 
the origin of such behavior in a more analytic way. 

Asymptotic behavior of one-electron Green's function- 
The correlation function ge{x,t) cannot be evaluated 
in a closed form, so that an analytic form of asymp- 
totic behavior will be of great help in understanding 
the nonmonotonic dependence of conductance on bias. 
First of all we note that the equal time Green's function 
ge{x,t = 0) is independent of interactions. However, we 
are mostly interested in the long time limit, so that we 
focus on the domain |a:/uct| <C 1. We write ge{x,t) in 
the following form: 



ge{x, t) — exp — const. + I{x, t) 
1 r°° dqe-^'"' 



I{x,t) 



q 



(28) 
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where const. = ^ 



dqe 



is a (infinite) constant, and 



(Pq{x,t) = qx + ujgt 



(29) 



is a phase function. In fact, the integral I{x,t) diverges 
at g = 0, and the divergence is cancelled by the above 
constant. In spite of this divergence, the form of Eq. ([28|) 
is more preferable for the asymptotic analysis. To avoid 
the divergence at q = 0, we employ the technique of 
dimensional regularizationifi: replace the factor 1/q of 
Eq.dlHl) by with a > 0, and take the limit a 

limit and extract the finite contribution. An infinity 
which appears in extraction is cancelled by the infinite 
constant mentioned above. 

In the long time /distance limit, the phase ip{x,t) be- 
comes very large in the generic domain, then we can use 
the stationary phase approximation.— In general, there 
exist two contributions to the asymptotic behavior of the 
integral of the type of /(a;, t): one from the ends of inte- 
gral interval (g = and q = oo, clearly the contribution 
from g = oo is negligible) and the other from the station- 
ary point(s) where dipq{x,t)/dq = 0. 

The contribution from the end point at q = can 
be obtained by integration by parts method.— Define a 
variable u as follows: 



u{q) 



qvc ln( 



' qa' 



(30) 



where -q = ^ef'-^/*-'^'''''/"'^. Changing the integration vari- 
able from q to u, the relevant integral becomes 



du\ 



-aq 



dq^ 

du [q{u)]^' 



(31) 



Let us consider the case where u > 0. Within loga- 
rithmic accuracy we have q{u) ~ ^^P'^ '71^]"^ ^^'^ ^ ~ 
-j- , / y with n* = e[(^/*)+^"l/^^-i. Performing the 

partial integration along imaginary axis^, taking a — > 
limit, and extracting the finite part, we obtain 

V t 

Icnd{x,t) -ln[— \n{riVct/a)] - z sign(i) 7r/2 (32) 
a 

The effect of LRCI within each channel is reflected in the 
double logarithmic correction. ln[ln{riVct/a)]. Note that 
this correction also depends on position li, I2 through 77. 
In terms of the Green function ge{x,t), the end point 
contribution is roughly l/[thit]^^'^ , which evidently can- 
not cause the nonmonotonic behavior of flux dependent 
conductance ct$ as shown in Ref. Q . 

Next we turn to the contribution from stationary point. 
The condition for the stationary point is 



dipq{x, t) 



dq 



^ ^ t 



dujq 

dq 



,\n^, (33) 
qcQ 



where ^ = Ce^""/"'^-'"^, and qc is determined by x/t. For 
the condition Eg. ([55]) to be satisfied for very long time 
limit (namely |a;/wci| <C 1) there should be a point where 



^^\q=q^ = 0. In other words, the frequency should at- 
tain a local maximum or mininum at finite momentum. 
Therefore if tOq is a monotonic function of q, then there 
will be no stationary point, so that the nonmonotonic be- 
havior of visibility would not appear. Another dispersion 
^q = voq — hq^ with h > which also shows nonmono- 
tous behavior satisfies the local maximum condition. It 
is easily seen '^gi^'*"' \q=qc ~ ~lt Then the standard 
stationary phase approximatioi>iS gives 



Ist&{x,t) 



27r 



\t\qcVc 



-if sgnt -i¥',(a;,t)|,=<j^ 



(34) 



In the long time limit |/cnd(a;,^)| ^ |^sta(a;, t)|. Were it 
not for the \/q singularity a,t q — Q, the stationary point 
contribution {l/\fi) would dominate the singularity- free 
end point contribution (1/i). Ea. ([M)l explains the os- 
cillating behavior of Green function at long time tail as 
found in Fig. 4 of Ref.fTi]. Now the asymptotic form of 
the correlation function ge{x,t) is given by 



9e 



ln(77z;,t/a)]i/^ 



(35) 



Discussions- The origin of the nonmonotonic behav- 
ior of visibility found in Ref. 7] is the oscillating tail of 
one-particle eletron Green function. In this report we 
have found that the oscillating tail is due to the specific 
momentum dependence of collective excitation. The in- 
vestigation of the detailed form of the visibility requires 
the numerical integration^, but the essential features can 
be understood qualitatively using Eg. ([55]) . Taking the 
typical time as t ~ l/\eV — i?c|, (A = hvc/a) one can 
estimate the integral of Eq. ((22| to obtain 



cr<i) 



\eV-EJ\n 



\eV-EJ 



(2/z.)-2 



' A ' 



(36) 



where (fidV) is basically the phase function ipq evalu- 
ated at the stationary point, and it slowly varies as a 
function of bias. The factor cos{(pc{V)) is responsible for 
the nonmonotonic behavior. At integer filling v = 1 the 
presence of energy scale E^. is not so pronounced because 
the exponent of prefactor of Eg. ([5S)) vanishes. Compared 
to Ref.f?!, the most dramatic difference with QH case 
(j/ < 1, for electron tunneling) would be the strong sup- 
pression of amplitude near eV ^ Ec- As for the quasi- 
particle tunneling (at v < 1), the time integral of Eg. p^ 
diverges in the long time limit, which implies that MZI is 
entering strong coupling regime at low temperature. It is 
well-known that the strong coupling regime of QH point 
contact is well-described by electron tunneling pictureJ^ 
At high temperature, the divergence in long time limt 
is cut by 1/T. Thus as temperature decreases, we can 
expect a crossover from quasiparticle tunneling to elec- 
tron tunneling, and at the same time the interchannel 
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Coulomb energy scale changes from ly^Ec to Ec- All these 
features can be readily checked experimentally. 

We have to note that the results of both Ref.Q and 
the present report have been obtained in the second order 
perturbation with respect to tunneling amplitudes. As 
such both results do not seem to explain the experimental 
results of the periodic behavior of visibility. However, 
the present report indicates the existence of a series of 
energy scales nEc {n = 1,2,3,.. .), and it clearly suggests 
that we need the higher order perturbations to reveal a 
possible underlying periodic structure of visibility. 
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